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Abstract 

The super Yang-Mills duals of open strings attached to maximal giant gravitons are 
studied in perturbation theory. It is shown that non-BPS baryonic excitations of the 
gauge theory can be studied within the paradigm of open quantum spin chains even 
beyond the leading order in perturbation theory. The open spin chain describing the 
two loop mixing of non-BPS giant gravitons charged under an su{2) of the so(6) R sym- 
metry group is explicitly constructed. It is also shown that although the corresponding 
open spin chain is integrable at the one loop order, there is a potential breakdown of 
integrability at two and higher loops. The study of integrability is performed using 
coordinate Bethe ansatz techniques. 



1 Introduction and Summary 



In the present paper we study the mixing of non-BPS giant gravitons in A/'=4 SYM beyond 
the leading order in the 't Hooft couphng. We pay special attention to open-string like 
operators coupled to maximal giant gravitons built out of two of the three complex scalars 
Z, W at two loop order. We show by an explicit computation that the two loop mixing of 
such operators can be described within the paradigm of open spin chains and we construct 
the two loop open spin chain. We thus extend the work of Berenstein and Vazquez |31 E] 
who had shown that the one loop mixing of non-BPS giant gravitons can be described by 
integrable open spin chains. Our analysis suggests that although the spin chain paragigm 
holds beyond the leading order in perturbation theory, the two loop open spin chain is not 
integrable. 

The discovery of integrability in A/'=4 SYMpLj has lead to various remarkable insights in 
the present day understanding of A/'=4 SYM and many of its less and non-supersymmetric 
cousins. The fact that in 't Hooft's large limit the mixing of single trace operators of the 
gauge theory can be studied and the spectrum of anomalous dimensions of the gauge theory 
computed by diagonalizing an associated one dimensional quantum spin chain, which, to the 
extent that it has been possible to check by explicit computations, appears to be integrable 
opens up the possibility of detailed higher loop studies of the gauge theory. Integrability in 
superconformal Yang-Mills theory was first discovered at the one loop level in the seminal 
paper of Minahan and ZarembopP who showed that the mixing of single trace operators 
constructed out of the six (real) scalars of the gauge theory is described by an integrable 
S0{6) spin chain. This result has subsequently been generalized and supersymmetrized to 
incorporate the mixing of single trace operators involving any of the SYM fields at the one 
loop level|2]. At this level, the mixing is described by a supersymmetric closed spin chain, 
which nonetheless is integrable. 

The status of integrability, beyond the one loop order is also very promising. Higher loop 
integrability in the gauge theory, which was first glimpsed in[^ has since been understood 
in much grater detail. It was shown by StaudacherjSj that higher loop integrability; at 
least in some closed sub-sectors of operator mixing such as the su{2) sector built out of two 
complex scalars, the sl{2) sector involving a complex scalar and a single covariant derivative 
and the S'u(lll) sector involving a complex scalar and a fermion, can be made manifest 
as one can write down a set of Bethe equations whose solutions determine the spectra of 
anomalous dimensions of the gauge theory and also produce the factorized scattering matrix 
of the closed spin chains that describe these sectors. The insights gained from the Bethe 
equations that determine the spectra in these various closed 'two-component' sectors have 
been utilized by Beisert and Staudacher in|nj where a remarkable all loop Bethe ansatz 
was proposed. Assuming the correctness of this Bethe ansatz, one is lead to the extremely 
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pleasant situation where the all-loop large spectrum of anomalous dimensions of all single 
trace operators of = 4 SYM is described by an underlying integrable quantum spin chain. 
Though, one does not yet know the Hamiltonian of this complete supersymmetric spin chain 
one can write down the all loop Hamiltonian in some less supersymmetric sectors. For 
example in the su{2) sector, an all loop Hamiltonian that reproduces the proposed su{2) 
Bethe ansatz and the appropriate dispersion relations can be written down in terms of a 
strong coupling expansion of the Hubbard model[H|. Thus, as far as the large dynamics 
of single trace operators of the gauge theory is concerned, one has many indications that it 
is governed by an integrable (closed) spin chain. 

Apart from single trace operators, the gauge theory also contains baryonic operators, 
whose engineering dimensions are of order A^. These operators are constructed out of de- 
terminants and sub-determinants rather than tracesfSl UHl E|- Such baryonic operators 
are of particular interest in the context of the AdS/CFT correspondence as they provide 
one with the gauge theory dual of non-perturbative D brane excitations i.e the so called 
giant gravitons of the dual string theorv [TH| IT^ OH 12m I22j . The simplest example of a giant 
graviton is provided by a determinant operator constructed out of one of the three complex 
scalar fields Z 

This operator being built out of only one of the complex scalars is a BPS operator and does 
not undergo renormalization. The gravitational interpretation of this operator is that of a 
maximal D3-brane that saturates the angular momentum bound [.I^ 1^1^ 1^ 1^. Non-BPS 
excitations in the background of a giant graviton would correspond to open string excitations 
in the presence of the brane. The gauge theory provides us with a simple enough prescription 
for constructing the boundary theory dual of the open string degrees of freedom. The open 
string ending on the D3 brane can be built by replacing one of the Zs in the determinant 
with a matrix product of local super- Yang- Mills fields [211 A typical non-BPS su{2) 
operator would look like 

e^'-'''ej,...j^Zi^ ■ ■ ■ Zi^zl {WWZZZWZWZ ■ ■ ■ ZZWZWy.^ . (2) 

Such operators, by virtue of not being protected by any non-renormalization theorems, would 
undergo operator mixing and hence have non-zero anomalous dimensions. The spectrum of 
anomalous dimensions of such operators gives one the gauge theory dual of the open string 
spectrum. Computing the spectrum of anomalous dimensions of such non-BPS maximal 
giants is the problem that we shall concern ourselves with in the present paper. As with 
various other issues concerning the AdS/CFT correspondence, it is important to analyze the 
problem in the large A^ limit before addressing the problem in all its generality. 

As it stands, it is not at all obvious that the computation of anomalous dimensions of 
giant graviton like operators simplifies in any way in the large A^ limit. Processes that naively 
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appear to be sub-leading order in may in fact turn out to be of order one as various ^ 
suppressions may be compensated for by the 0{N) size of the operators. The fact that 
such counter-intuitive turn of events can come about in the study of Yang-Mills theories 
is rather well known from the study of the A/'=4 SYM in the plane wave limit. There too 
one studies long operators that have large bare conformal dimensions and R-charge, J, and 
one has the result that if the ratio is held fixed then the usual notions of planarity do 
not apply any more 7j. Non-planar contributions turn out to be of 0(1) in that limit. The 
study of operator mixing can still be reduced to a quantum mechanical problem, though the 
underlying quantum mechanics that one is faced with is not that of a quantum spin chain 
but that of a non-local many-body problem. Thus one of the pleasant features of the gauge 
theory, namely integrability is no longer present in this novel-large N limit^. 

Though we shall not aspire to construct any new large limit where some parameter 
(other than A^, such as the R-charge) is allowed to approach a critical value, and work 
completely within the framework of 't Hooft's large limit, the baryonic operators that we 
shall work with are indeed of order A^. Thus we shall keep the lessons learnt from the study 
of the plane wave limit as a valuable caveat. From the studies of the dynamics of non-BPS 
giant gravitons of the gauge theory performed so far, the status of integrability is ambiguous. 
It was first shown in|21] that, at least at the one loop level, the anomalous dimensions of 
non-BPS giants are not of order one i.e although these operators have large bare conformal 
dimensions, their anomalous dimensions are small. Subsequently, more detailed one-loop 
analysis of non-BPS operators involving all six (real) SYM scalars was carried out ini4j. In 
that paper, one loop mixing of operators such as 

where ^*s can be any of the six scalars, was carried out. It was shown that the computation of 
one-loop anomalous dimensions of such operators can be reduced to the problem of diagonal- 
izing the Minahan-Zarembo so(6) quantum spin chain, but with open boundary conditions. 
The open boundary conditions emerge naturally in the problem as the operators of interest 
are not traces, and the interaction of the scalar fields that form the chain(the matrix product 
of the (p fields in the above operator) with the Z fields forming the brane lead to the open 
boundary conditions. The outcome of the analysis performed in^^ was the discovery that the 
open boundary conditions do not spoil the integrability of the Minahan-Zarembo spin chain. 
It was shown in that paper that the bulk and boundary S matrices for the relevant so{6) 
spin chain satisfy the Yang-Baxter and reflection algebras: thus a necessary condition for 

^In various other scenarios as well, where one takes into account non-planar splitting and joining of strings 
the notion of integrability, if any, appears to be far from obvious. We shall refer to |13[ ITl) for some recent 
developments in this direction. 
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integrability is indeed satisfied at the one-loop level. However, already at the one loop level, 
the situation is different if one considers non-BPS non-maximal giants, which correspond to 
operators such as 



with A; > 1. w stands for the matrix product formed out of two complex SYM scalars 
i.e. it is the su{2) spin chain. It was shown injB] that the mixing of non-maximal giant 
gravitons of the above type formed out of two complex scalars can be formulated in the 
language of quantum spin chains, at least at the one loop level. The spin chain is nothing 
but the celebrated Heisenberg spin half model, but the boundary conditions generated by 
non- maximal giants are of a novel sort. The boundary interactions result in the spin chain 
being a dynamical one, i.e, the boundary interactions lead to fluctuations in the length of 
the chain. Though, it was possible to find the exact ground state of the spin chain in|3] 
whether or not the chain is integrable is still an open issue^. Thus, integrability in N=4 
SYM is not at all obvious once one ventures outside the realm of single trace operators. 

Given this status quo, it is probably a reasonable goal to try and examine whether the 
spin chain interpretation of mixing of non-BPS giant gravitons holds beyond the one loop 
level and if it does it is of course interesting to investigate whether or not they can be 
diagonalized by Bethe ansatz techniques that have been so successful in the study of higher 
loop integrability of single trace/closed spin chain operators. In the present paper we shall 
take a step in that direction. The organization of the paper and the principal results reported 
in various sections are as follows. 

In the next section we shall establish that the dynamics of non-BPS excitations in the 
background of a maximal giant graviton can be studied within the paradigm of open spin 
chains at the two loop level. In other words the mixing of non-BPS maximal giants will 
be shown to be described by open spin chains at the second order in the 't Hooft coupling. 
The analysis will be performed in two steps. We shall first resolve various combinatorial 
issues to work out the proper normalizations for various determinant and sub-determinant 
like operators. The guiding principle used in this exercise is that the various operators are 
required to have a free field two point function of order one. Once the combinatorial issues 
are settled, we shall proceed to study the mixing of these non-BPS giant gravitons. To that 
end we shall use the results of Beisert et al[3^ where the two loop dilatation operator for the 
su{2) sector of the gauge theory was computed at finite A^. Isolating the 0(1) part of the 
result of the action of the finite dilatation operator on the properly normalized operators 

^The dynamical nature of the spin chain does not by itself rule out integrability. Even in the sector 
of closed chains e.g in the closed su{2\3) sector dynamical chains do appear [TT] and they appear to be 
integrable|12j. However, it is not known if the dynamical spin chains of interest in the sector of giant 
gravitons of the gauge theory enjoy similar integrable features. 




(4) 
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gives us the details of mixing of non-BPS giants at large A^. We shall be able to show that 
at the two loop level non-BPS giants of the type© are closed under action of the dilatation 
operator. We shall also be able to interpret the process of mixing of such operators within 
the language of open spin chains; the open spin chain being the matrix product of various 
su{2) SYM fields that are contracted with the last pair of indices of the epsilon tensors inQ. 
We shall also be able to show that no new boundary conditions are generated at the two 
loop level. Thus the boundary conditions are the same as the ones discovered by Berenstein 
and collaborators in^, which amount to the fact that the first and the last fields in the 
open chain© cannot be Zs. These boundary conditions emerge out of the realization that 
if a Z field is placed at either end of the open spin chain, then the operator factorizes into 
a determinant and a trace, which is lower order in This fact may most easily be seen by 
going to a basis in which Z is diagonal. Thus the boundary conditions for the problem at 
hand are not exactly the so called free boundary conditions, which would have corresponded 
to the case of the first and last spins being allowed to take any one of the two allowed values, 
but they are very similar. 

Having a spin chain interpretation of the mixing of non-BPS giants at two loops is a 
gratifying outcome. The one loop open spin chain was already shown to be integrable inj4j, 
and our analysis suggests that the study two loop operator mixing can be regarded as the 
analysis of a perturbation to an integrable open spin chain. The obvious question that 
comes to mind is whether or not the two loop spin chain Hamiltonian can be diagonalized 
by an appropriate Bethe ansarz. That, unfortunately does not appear to be the case. Before 
indulging in further details, it is probably worth emphasizing that even in the absence of 
integrability in the sense of the absence of a Bethe ansatz and a factorized S matrix, the 
fact that the two loop operator mixing can be described by open spin chains is a major 
simplification in the study of the open string/ giant graviton sector of the gauge theory. That 
is so because, as far as the computation of the spectrum of anomalous dimensions of operators 
that correspond to spin chains of finite length is concerned, the spin chain interpretation 
always reduces the problem to the process of diagonalization of a finite dimensional matrix, 
which can be achieved numerically. However, the lack of a two loop Bethe ansatz is a major 
drawback if one intends to study long coherent open string states constructed from the gauge 
theory and compare with the world-sheet analysis as was done in^|. 

Having established the open spin chain description of the mixing of non-BPS giants, 
we shall devote the rest of the paper to a careful analysis of the integrability or the lack 
thereof of the two loop open spin chain. Although integrability of a given spin chain, even 
in a perturbative sense, is in general quite hard to establish, we shall take the pragmatic 
approach of StaudacherjH] and declare the spin chain to be integrable if one can successfully 
write down a set of Bethe equations for the spin chain. In the approach pioneered in0, the 
basic idea was to study the two magnon scattering problem, and find the appropriate two 
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body 5* matrix. Assuming factorized scattering, the two solution of the two body problem 
was then used to find the exact (in a perturbative sense) eigenstates, i.e the Bethe eigenstates 
of the general many body problem. An explicit construction of the scattering eigenstates of 
the su{2) dilatation operator and the corresponding Bethe equations, accurate to the third 
order in the 't Hooft coupling was also worked out by Staudacher in the same paper Bj- 
As far as the closed chains of the su{2) sector is concerned, the main complexity lay in 
the construction of the two magnon scattering matrix and the Bethe equations describing 
the two magnon scattering. As long as there is no global obstruction to the construction 
of the two particle S matrix, one could assume factorized scattering and lift the solution 
of the two body problem to the case of an arbitrary number of magnons. One could then 
carry out further checks to ascertain that the solutions of the general many body problem 
do in fact reproduce the correct spectrum. In the su{2) sector, for instance, the spectrum 
obtained from the Bethe ansatz|B^ does match the spectrum obtained from the numerical 
lattice simulations of the corresponding spin chain lOj. In a sense, the notion of integrability 
employed in this approach is the factorized nature of the scattering matrices for three or 
more magnon scattering processes. 

However, in the presence of boundaries, there is a second notion of three body scattering, 
which is present even in the analysis of the two magnon problem. The presence of the 
boundary can lead to three body processes even if one is dealing with the dynamics of two 
magnons. Integrability in the presence of boundaries hence implies that these three body 
interactions take place in a way that is compatible with the separation of the scattering 
matrix into a bulk and a boundary part. Thus, even from the point of view of Staudacher's 
PABA apprach, the criterion of integrability is more stringent in the presence of non-periodic 
boundary conditions. We shall be able to demonstrate, in an explicit fashion that such three 
body processes, though absent at the one loop level, do start contributing at O(A^). The 
algorithm that we employ to do this is the following one. We compute the bulk and boundary 
S matrices from the study of the asymptotics of the two magnon problem. In other words, 
the study of the magnons scattering off each other far away from the boundaries gives us 
the bulk two body S matrix of the spin chain and the situation when one of the magnons 
scatters off the boundary (while being far enough away from the other magnon) gives us the 
boundary S matrix. From the knowledge of these scattering matrices, we can write down 
a set of Bethe equations. We then proceed to show that these Bethe equations follow from 
the requirement for a general ansatz for the two particle wave function to be an eigenstate 
of the spin chain Hamiltonian. Finally we demonstrate explicitly that the two body wave 
function that leads to the Bethe equations is not an eigenstate of the Hamiltonian. We 
interpret this to be evidence in favor of a breakdown of higher loop integrability in the open 
string sector of the gauge theory. Though the Bethe equations do not appear to be solutions 
of the two body Schrodinger equations, we can nevertheless utilize them to construct some 
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exact eigenfunctions of the Hamiltonian. This construction is discussed at the very end of 
the third section. 

From the point of view of the spin chain, we can also understand the potential breakdown 
of integrability in another way. The two loop dilatation operator that describes the mixing 
of non-BPS giants is the same as the one for closed chains as long as one stays away from 
the boundaries, i.e the bulk dilatation operator can be embedded in the Inozemtsev spin 
chain^Q. To the best of our knowledge, the only boundary conditions, other than the 
periodic ones, that are compatible with the integrability of the Inozemtsev chain are provided 
by reflective boundaries (see for example [23 l2S»i22j). The boundary conditions that emerge 
out of the gauge theory analysis are not of that type, and as mentioned earlier, they are 
similar in spirit to the free boundary conditions. The Inozemtsev chain is not known to be 
integrable in the presence of free boundary conditions, and hence it is perhaps not surprising 
from this point of view that the two loop dilatation operator in question in the present paper 
does not appear to admit a consistent Bethe ansatz. 



2 Various Operators: 

The generic class of operators of interest to us are the ones that describe non-BPS excitations 
in the background of a maximal giant graviton. As mentioned earlier in the introduction, 
maximal BPS giants are protected operators formed out of determinants of a single complex 
scalar Z, and they define the vacuum state for the sector of interest to us. 

\0>=e^:.ZZ^...Z;- (5) 

It is understood that 

Jl-jN _ . Jl-jN ((W 

A typical open string excitation corresponds to a state such as 

IZ-- '7™-,,, 1,, ^ il---ikik+l---iN yjk+l 7^>' + ^ ( nn \ik+m + l ( \jN ( J\ 

\k,Z ,Wi,---WN~k~m>- ei^...idk+l-jN^ik+, • ■ • (W'lJife+^+i ■ ■ ■ l^TV-fc-mji^ 

The WiS are matrix products of local Yang-Mills fields i.e they are the open string degrees 
of freedom. In the above operator, the indices of the epsilon tensor {k of them) which are 
contracted with each other indicate the departure from the maximal R caharge bound, while, 
replacing a certain number of the Z fields with the open strings, Wi, signal the departure form 
the BPS bound. As has been argued and demonstrated in[25j, operators of the above kind 
do satisfy orthonormality conditions and it is possible to build a large Hilbert space using 
these open string degrees of freedom once the states are properly normalized. Rather than 
review the detailed combinatorial work presented inj2S], which deals with the generic class 
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of operators to which the above example belongs we shall focus on the particular operators 
of interest to us. We shall be interested in operators of the type 

\Z''~\w> (8) 

which gives us the gauge theory description of a single open string ending on a maximal 
giant graviton. 

w] = {wiW2 ■ ■ ■ Wl)) ■ (9) 
Each of the 'bits' of the open string Wi can take on two values Z and W . 



2.1 Normalization: 

To be able to develop a consistent large expansion one needs to normalize the operators 
in a way that lends itself to a well defined large limit. We shall chose to normalize 
the operators by the free-field two point functions. It may be seen by a straightforward 
computation that 

(z^-i(x),^(x)|Z^-i(0),MO))~ [(iV-l)!]3iV^+i^^^ (10) 

Hence, it is reasonable to normalize the states as 

, ^ \Z'^-\w> (11) 

y^[(iV- l)!]3A^^+i 

As mentioned in the introduction, with this normalization one is led immediately to dirichlet 
boundary conditions for the open spin chain w. It may be easily verified that if either wi 
or Wl is equal to Z, then the state factorizes into a product of a determinant and a single 
trace jH 123- With our chosen normalization, such a state would turn out to be sub-leading 
order in the large A^ limit. Thus, one has the Dirichlet boundary conditions that, neither 
the first or the last impurity fields can be of the Z type. 

2.2 The Open Spin Chain: 

The two loop dilatation operator of A/'=4 SYM in the su{2) sector is[Sj 

D = D, + ^D, + ^^D, + 0{)<^) (12) 

where 

Do = tr [A^^Az + A^'^Aw) (13) 
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tr[A^',A^'^][Az,Aw] 



(14) 



D2 



tri-[Az,Aw][A^MAz,[A^^A^'^]]] 

+^[Az,A^][A^'^,[AwAA^z,AU]] 
+ [A^^,A^'^][Az,Aw]), 



(15) 



and 



A 



9ymN 
87r2 



(16) 



It is understood, that one has mapped the local composite {su{2)) operators to states of the 
matrix model Hamiltonian. For example, 



and that operator mixing is mapped to the action of the dilatation operator, thought of as a 
matrix model Hamiltonian, on various states like the one depicted above. When the dilata- 
tion operator acts a giant graviton like state such as \Z^~^,w > the following possibilities 
may occur. 

(1:) All the annihilation operators of Di and D2 may act on w. 

(2:) Some of the annihilation operators may act on w and others on the Z fields that form 
the brane. 

In the first instance, one has the standard result that the action of the annihilation opera- 
tors on consecutive fields present in w can be interpreted as the action of a quantum spin 
chain Hamiltonian on a spin chain state jHl ED]- The action of the annihilation operators on 
non- consecutive fields can be shown to be sub- leading order in Such actions always spit 
up the spin chain into multiple chains, thus transforming the state with a single string exci- 
tation in the D3 background to one involving multiple open string excitations. Such states 
are mutually orthogonal (see[2S] for the combinatorial details). Hence they can be ignored 
for the purposes of large N computations. The spin chain that results from the 'bulk' action 
of the dilatation operators on consecutive fields present in w is the same as the one that one 
has in the study of single trace operators [3 , only the boundary conditions are different. As 
mentioned previously, the boundary conditions consist of the condition that neither the first 
or the last fields of the spin chain can be a Z. 

When the dilatation operator acts on two Z fields from the brane, the result can be easily 
shown to be sub- leading order in j^. However, special care is required to analyze the result 
of the action of the dilatation operator on two consecutive fields in w and a single Z field 



Tr{ZZWW---W) 



rr(At^At^Atw^...At^)|o>, 



(17) 
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from the brane. A careful analysis shows that there are no 0(1) contributions from such 
actions of the dilatation operator. Indeed, all the terms so generated are given below along 
with the leading powers of jj: by which they are suppressed. In the following list, wi stands 
for a single field, W2 represents a generic 'word' formed out of two fields and Wc is a word of 
length — 2: \w\ being the length of w. The second set of states multiplying the fifth and 
lower terms in the table are single trace states which result upon the action of the dilatation 
operator on the giant graviton state. 

Wi, WcW2 > 

\Z''-\ws,w,> 

\Z^~'^, W2, WiWc > 

\1,Z'^-\W,> \W3 > 
\1, Z'^~^ ,W2Wc > \Wi > 
\1, Z'^~^ ,WiWc > \W2 > 

\Z'^~'^,W2,Wc > \Wi > 
\1, Z'^~'^ ,W2Wc > \Wi > 

\Z'^''^,Wi,Wc > \W2 > 

We can now write down the formula for the two loop large open spin chain. As mentioned 
in the preceding discussion the two loop Hamiltonian for the open string sector is the same 
as the one for single trace operators, the only thing that is is different is the boundary 
condition. To incorporate the boundary conditions in the formula for the Hamiltonian we 
shall use the operators Qi{Z), which were introduced in|l]. These operators are defined by 
their actions on the spin chain states as foillows: 

Qf\---Z---> = 

Qf\...W---> = \---W--->. (18) 

In the above equations, Z,W are located at the 'i'th lattice site. i.e. the operator Qi{Z) 
checks if the 'i'th spin is a Z. If it is so then it annihilates the entire state. Using these 
operators, the open spin chian Hamiltonian can be written as 

H = J2i^-4g)QfQj{I-Pi,i+,)QfQ^j 
1=1 

+9j2QfQ5ii-PiH2)QfQ5 

i-i 



1 

1 

J_ 

N 
1 

Vn 
1 
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To get to the above form, we have divided out the formula for the Hamiltonian by A and 
used 

9=\. (19) 

It is of course assumed that the length of the chain is J and the Hamiltonian acts on all the 
spin sites. Since the first and the last sites cannot have Z fields, one could write down the 
Hamiltonian is a slightly different form as well. We could take a state of the generic form 

\W---W>, (20) 

where one could have impurity {Z) fields inserted between two W fields. In the previous 
formulation, the spin chain included the boundary fields, but we could just as well modify 
the definition of the spin chain to include only the fields contained inside the two boundary 
W fields, and write the effect of the interaction of the boundary fields with the spin chain 
separately. In other words, we shall define the second field of w to be first site of the spin 
chain and the penultimate field of w to be the last site. Thought of in this way, the spin 
chain Hamiltonian reads as 

H = X: (1 - 4^)(/ - Pm+i) + 9 eV - ^M+2) 
1=1 1=1 

(1 - Ag)qf + gql + (1 - Ag)ql + gql_, (21) 

It is to be understood that 

gf = (l-Qf), (22) 

i.e qf checks if the field at the ith site is Z, and if it is so, it acts as the identity operator. 
In case this condition is not satisfied, qf annihilates the spin chain. L is now to be regarded 
as the effective length of the spin chain i.e the number for fields between the boundary W 
fields i.e. 

L=\w\-2. (23) 



2.3 Bethe Ansatz: 

Let us begin by deriving the one loop Bethe ansatz for the open spin chain. Short ranged 
spin chains of the Heisenberg type with various boundary conditions have been solved using 
both the coordinate and algebraic Bethe ansatz techniques in the past jHH IS21 iHB] ■ In the 
present paper, we shall adhere to the coordinate Bethe ansatz approach as it generalizes 
easily to higher loops. An excellent discussion of coordinate Bethe ansatz solutions to short 
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ranged spin chains can be found in[^ |32|, and we shall adhere to the approach presented 
in these papers in this present analysis. 

The one magnon state |\E'i > can be taken to be of the form: 

>= > . (24) 

If one takes the state to be a superposition of incoming and outgoing plane waves, with 

^(x) = Aiky^"" - y4(-fc)e-*'=^, (25) 
then it is quite easy to see that the state is an eigenstate of the Hamiltonian with energy 

E(A;) =4sin2 . (26) 

It is also straightforward to verify that for the plane wave state to be an eigenstate, the 
following boundary conditions, which are implied by scattering at the boundary, must hold. 

$(0) = 

^(L + 1) = 0. (27) 
These conditions can also be written down as 

a{-U)A(}i) - a{}z)A{-U) = 

P{k)A{k) - P{-k)A{-k) = 0. (28) 
The above conditions lead to the following compatibility condition 

a{k)(3{k) = a{-k)(3{-k) (29) 

with a, (3 given by 

a{-k) = 1 

(3{k) = e^(^+l)^ (30) 

The solution to the compatibility condition sets the allowed range of values for the momenta 
to be 

k=-^. (31) 
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This is the complete solution to the one magnon scattering problem, which gives us the 
one-magnon spectrum ()26l31|) and the boundary scattering matrix ^^^^ . 

Having determined the boundary S matrix by solving the one-magnon problem, one can 
now proceed to study the two magnon scattering problem. The two magnon state can be 
taken to be of the form 

1^2 >= > (32) 

x<y 

with 

■^ix,y) = Y.aip)Aih,k,)e'^'^^^-''''-'\ (33) 
p 

The above sum is over all permutations and negations (p) of the momenta. a{p) = ±1 and 
it changes sign when the signs of when one of the signs of the moments are changed or when 
the momenta are permuted in pairs. This two magnon state is an eigenstate of the spin 
chain Hamiltonian provided the following conditions are satisfied. 

A{ki, k2)s{ki, k2) - A{k2, ki)s{k2, ki) = 
A{ki, k2)a{~ki) - A{-ki, k2)a{ki) = 

A{ki, k2)(3{k2) - A(A;i, -k2)(3{-k2) = (34) 

In the first equation 

s{ki, k2) = l- 2e^^ + e'^^'+^^\ (35) 

The three equations above require some explanation. The first of the three equation comes 
from the requirement that two body scattering, that takes place when the impurities appear 
next to each other, happens in a way such that the state remains an eigenstate of the Hamil- 
tonian. The other two equations follow from scattering at the boundary, i.e by considering 
the scenario that one of the particles is at the boundary and by requiring the coefficient 
of this state to be such that it still remain an eigenstate of the Hamiltonian. In a sense 
we can see, directly from these equations that the boundary conditions are compatible with 
integrability at the one loop level. That follows from the observation that the boundary 
scattering equations involve only one of the two momenta. In other words, the boundary 
does not introduce any new 'three body' interactions. For example, if one considered the 
special case when the two particles sit next to each other at the boundary, there could in 
principle have been potential three body interactions; the third body being the boundary. 
However, the above conditions tell us that this is not the case at the one loop level, and the 
two boundary equations are all that one needs for all positions of the second impurity. We 
shall soon see that this is not the case at the two loop level. 
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The situation is more transparent in position space. For any pair of well separated lattice 
sites away from the boundary, 

H\x,y{> X + 1) >= 4:\x,y > —\x — l,y > —\x + l,y > —\x,y — 1 > —\x,y + 1 > (36) 

With the choice of the pane wave ansatz, and this action of the Hamiltonian, it is clear that 
the coefficient of \x,y{> x + 1) >, once the Hamiltonian has acted on it, is (4sin^ (J^^ + 

4sin^ (J^^)'^{x,y). However, for the coefficient of \x,y{= x + 1) > to come out right, one 
has to impose the condition 

2^(x, X + 1) - ^(x + 1, X + 1) - ^(x, x) = (37) 

which leads to the scattering condition; the first equation in Once this equation is 

satisfied, the effects of the boundary impose the further constraint 

^{0,y) = 

*(x,L + l) = (38) 

which translate into the second and third equations in ()34j] respectively. What is notable 
is that once the two particle scattering equations are satisfied, for some choice of \l/(x,?/), 
the boundary conditions on \E'(x,?/) are the same as the one magnon case. For instance, the 
left boundary condition (the first of the two equations above) is the same as the condition 
satisfied by ^(x) in the one- magnon problem for all values of y. Most importantly, the 
special situation when x = l,y = 2 does not introduce any new conditions on the wave 
function. Thus the bulk and boundary scattering problems separate out. This is nothing 
but an explicit (and perhaps pedagogical) way of realizing that the boundary conditions are 
compatible with integrability of the spin chain. 

Having the three equations at out disposal, we can now derive the Bethe equations that 
determine the momenta. As in the one magnon case, the Bethe equations are derived as the 
compatibility conditions implied by the three equations. Various ways of expressing A{ki, /C2) 
in terms of A{k2,ki) lead to the requirement |32j 

«(fci)/3(fci) ^ S{-k,,k,) 
a{-k,)(3{-k,) S{k,,k,) ' 

where 

S{k,,k2) = sik,,k2)s{k2,-ki) (40) 

There is a similar equation where ki is replaced by ^2 and vice-versa. These results can be 
generalized to the case of an arbitrary number of magnons. In the case of m magnons the 
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(39) 



Bethe equations read as 



n 



m 



S{ ki, kj) 



(41) 



a{-k,)f]{-k,) 



S{ki, kj) 



3 The Two Loop Analysis: 

Bethe ansatz for the su{2) open spin chain decribed in the previous section is essentially 
implied in the analysis carried out in[4r where the one loop bulk and boundary S matrices of 
the full so(6) spin chain were computed. The above analysis reconfirms and reproduces the 
su{2) sector of the results reported in|ljfrom a coordinate Bethe ansatz point of view and 
also sets up the techniques and conventions that will be applied to the two loop analysis, to 
which we shall devote the rest of the paper. The two loop large dilatation operator is 



As was the case in the analysis of single trace operators, we do not expect the two loop spin 
chain to be exactly integrable. However, in the analysis of closed spin chains, it was possible 
to establish a perturbative integrablity of the large dilatation operator. In other words, 
it was possible to construct a Bethe ansatz for the eigenstates of the two loop Hamiltonian 
which was accurate to 0{g)^. It is thus reasonable to examine if the two loop large A^ 
dilatation operator exhibits this perturbative notion of integrability in the background of 
giant gravitons. 

The method of the perturbative assymptotic Bethe ansatz (PABA) advocated by Stau- 
dacher injH] provides one with a very direct and transparent method for checking whether 
or not a given spin chain, whose Hamiltonian can be expressed as a power series in some 
coupling constant (in our case g), admits of a Bethe ansatz order by order in perturbation 
theory. This method was successfully applied to the study of the mixing of single trace su{2) 
operators of the gauge theory up to three loops in perturbation theory. In the sector of single 
trace operators, the large A^ dilatation operator maps to closed spin chains with periodic 
boundary conditions. For such spin chains, there are no boundaries, and the magnons only 
scatter off each other. Establishing the existence of a PABA for the closed spin chains thus 
consists of a two step process. 

( 1:) Finding the two body S matrix for two magnon scattering while allowing for the 5" 
matrix to have an appropriate g dependence such that the two magnon Bethe state is an 
eigenstate of the Hamiltonian in the perturbative sense mentioned above. 




H=Y.{l-Ag){I- + ^7 E - ^M+2) 



1=1 1=1 



(1 - Ag)qf + gql + (1 - Ag)ql + gql_ 



(42) 
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( 2:) Generalizing the two body solution to a solution of the many body problem using the 
kind of reasoning used at the end of the previous section, which basically amounts to assum- 
ing factorized scattering for the problem. 

For the problem at hand, the only difference between the two loop Hamiltonian ()42|l and the 
corresponding spin chains describing the two loop dilatation operator for single trace states 
lies in the boundary conditions and the boundary terms(the second line in ()42|) ). Thus, 
when the magnons are far away from the boundaries, the scattering between the magnons is 
described by the same S matrix that appears in the case of the closed spin chains, which was 
already found in||Bj. Thus, for the purposes of the present analysis, we can borrow fromjH] 
the formula for the two loop bulk two magnon scattering matrix. However, for the present 
problem, we also need to perform a PABA to account for the scattering of the magnons from 
the boundaries i.e compute the boundary scattering matrix accurate to two loop order. As 
in the case of the one-loop analysis performed earlier in the paper, the boundary scattering 
matrix, will be found by solving the one magnon problem. To clinch the issue of perturbative 
integrabilty in analysis of the open spin chain would require a demonstration that the Bethe 
state formed by merging the two asymptotic solutions is an eigenstate of the Hamiltonian. 
To put it differently, we shall have to check that two magnon state that has as the bulk S 
matrix the one found in^, and as the boundary S matrix the appropriate solution of the one 
magnon problem is an eigenstate of the Hamiltonian. This unfortunately is not the case. We 
shall now proceed to substantiate these comments by explicit computations in the following 
sections. 

3.1 One Magnon Problem: 

As in the one-loop case, we shall first compute the boundary S matrix of the spin chain by 
solving the one-magnon problem. An important distinction from the one-loop analysis is that 
the Bethe wave function is not simply a superposition of plane waves. We shall also have to 
allow for some exponential decay of the wave function around the boundaries. The analysis 
below can be thought of as the two loop PABA for the one magnon problem. Although this 
is ostensibly a one-body problem, the presence of the boundary induces a two body effect 
through the interactions of the impurity with the boundary. Just like the two body problem 
in Staudacher's PABAjS] required the introduction of an exponential decay of the two body 
wave function for the Bethe ansatz to be consistent, the one magnon problem in the present 
analysis requires that the one body wave function be allowed to have an exponential decay 
around the boundaries. The Bethe ansatz that incorporates this requisite feature can be 
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summarized in the following form for the one magnon wave function: 

L 

1^ > (43) 



where 

- A{k)e'^'' - A{-k)e-'^'' + a{k)g'' + b{k)g^-'' (44) 

As is evident from the above formula, the third and the fourth terms account for the expo- 
nential decay of the wave function; the strength of the decay being given by the 't Hooft 
coupling constant. To avoid having to write factors involving the momenta, let us call 

Aiky^"" - A{-k)e-'^'' = i^{x) 
a{k) = 0(1) 

h{k) = 0(L) (45) 

As in the one loop analysis, when the impurity is sufficiently far away from the boundary i.e 
for any 2 < x < L — 2 we have the condition: 

(1 - Ag){2^p{x) - ijj{x - 1) - ip{x + 1)) + g{2ilj{x) - ij{x - 2) - il;{x + 2)) = E^{x) (46) 

with 

Eik) = 4sin2 - le^sin^ = E,ik) + gE^ik). (47) 

Special care is needed to analyze the situations x = 2, 1, L, L — 1. The coefficient of |x = 2 >, 
after the action of the Hamiltonian on the Bethe wave function is 

(2 - 4^)(2V'(2) - V'(l) - V'(3)) + g{2i;{2) - ^^(4)) - g<f>{l). (48) 

Thus we have the condition 

^(0) = 0(1) (49) 

Similarly, the coefficient of |x = 1 > is 

(1 - 4^)(2V;(1) - V(2)) + gm) - V^(3)) + 2#(1) (50) 
This leads to the condition: 

(1 - 4g)i;{0) + gi^i-l) - #(1) + 2#(1) = gE,{k)4>{l) (51) 
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where Ei{k) = 4sin^ ^|^. There are similar equations that follow from the scattering at the 
right boundary. These equations may be written in a way that is at par with the one loop 
equations (j28|l i.e 

A{k)a{-k) - A{-k)a{k) = 

A{k)(3{k) - A{-k)[3{-k) = 0. (52) 

determine the two loop left (right) boundary scattering amplitudes and are given by 

a{-k) = (l-2g) + ge-^'-ge''-gE,{k) 
P(k) = e'''^'^^'Uil-2g)+ge''-ge-''-gE,ik)). (53) 



As at the one loop level, the compatibility conditions a{k)(3{k) = a{—k)P{—k) determine 
the momenta to be 

k = (54) 

where 

N if 2gsm(k) 

9{k) = tan~^ ^ ^ 



l-2g)-gE,{k)\- ^^^^ 

It is understood that the expression for 6 is to be regarded as an approximate expression 
accurate to 0(5^). 0{k) gives us the two loop correction from the boundary scattering to the 
allowed values of momenta. This completes the Bethe ansatz for the one magnon probelm. 



3.2 Two Magnon Problem: 

Having solved the one magnon problem, we can now analyze the case of two magnons. Once 
again, we shall have to propose a Bethe ansatz, that accounts for the exponential decay of 
the wave functions when the magnons interact with each other or with the boundary. The 
appropriate PABA can be summarized in the following two body wave function 

1$ >= Y,'^ix,y)\x,y >= Y.ij{x,y)\x,y > +£/ ^ </)(!, j/)|l, y > +g ^ 4){x,L)\x,L> . 

x<y x<y y>l x<L 

(56) 

As in Saudacher's PABA[8J we shall let ip{ki,k2) be a superposition of plane waves, with 
some exponential decay that takes place when the magnons interact with each other, while 
(/), (j) will measure the amplitude of decays around the boundaries. More specifically, we shall 
have 

^(^1, ^2) = E a(p)e*('=i^i+^^"^) (A(A;i, k^) + B{ki, fcs)^?!"^-"^!) . (57) 

V 
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Once again, for two magnons at xi > 2 and L — 2 > X2 > xi + 2, the coefficient of \xi,X2 > 
after the action of the Hamiltonian is 

(1 - Ag){4:1p{Xi,X2) - i^iXi - 1,X2) - i^iXi + l,X2) - ^(Xi.Xs - 1) - ^{Xl,X2 + 1)) 

+g{AiJj{xi, X2) - ^{xi - 2, X2) - ip{xi + 2, X2) - ip{xi,X2 - 2) - ip{xi,X2 + 2)) 

= E{h,k2)i^{xuX2)= V (Asm'('^]-16gsm'(^]]ij{x,,X2) (58) 



At the special situation when X2 = xi + 2, one has the condition 

2'il){x, a; + 2) - + 2, x + 2) - ipix, x) = 0, (59) 
while when X2 = xi + 1, one has the condition 

gilpiXi + 2, X2) + i>{xi,X2 - 2) - ll){xi - 1, X2 - 1) - 1p{xi + 1, X2 + 1)) 

{Ag - l)(2V^(xi, X2) - V^(xi + 1, X2) - ^(xi, X2 - 1)) = (60) 

This part of the computation is the same as the PABA for periodic spin chains carried out by 
Staudacher in[Hj. The above conditions serve to fix the function B{ki, ^2) and to determine 
the bulk scattering matrix. However the boundary processes have to be treated separately 
which we shall now proceed to do. Now, for any X2, the condition obtained from looking at 
the coefficient of |2,X2 > is 

^(O,X2) = 0(1,X2). (61) 

One also has a similar condition at the right boundary; namely 

7/^(xi,L + l) = 0(xi,L). (62) 

Furthermore, for any X2 > 2 and Xi < L — 1 the coefficients of |1, X2 > and \xi, L > lead to 
following two boundary scattering equations respectively. 

A{ki, k2)a{-ki) - A{-ki,k2)a{ki) = 

A{ki, k2)P{k2) - A{k,, -k2)(3{-k2), (63) 

with a, (3 given bv fj^Hj) . Equations lUUl lUTl lUHj) completely fix all the undertermined 
functions in the problem. Indeed, they also lead to the Bethe equations for the momenta 

<^{h)m) -^s{-K,k,)S{k,,k;). (64) 



a{-ki)(3{-ki) 
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a and j3 are given by while the bulk S matrix is the solution of ()59|60|61|62|) . It is the 
same as the one presented injHl EI- 

Sfert = (65) 

with the phase shift (thought of as an expression accurate to order A for the present problem) 
being given by 

1 /p" 



*>(p)=2"t(f) 




1 + 8A sin^ \^\. (66) 



This situation is very much at par with what one had at one loop ()34p . The boundary 
interactions so far have not introduced any three-body interactions. The effect of the left 
(right) boundary scattering on '0(xi,a;2) produces the same algebraic equations satisfied by 
il){x) in the one magnon problem. Thus in the analysis so far, the boundary has not 'seen' 
the second particle. That is not surprising as we have not analyzed the only situation that 
can lead to a three body process, which is when both the impurities are near the boundary: 
i.e we need to study the coefficient of the state |1,2 > after the action of the Hamiltonian 
on the Bethe wavefunction. This is the only scenario (at this loop order) in which all the 
three 'particles', the two impurities and the boundary can interact with each other. Thus 
the spin chain will be perturbatively integrable, if the coefficient of |1,2 > does not lead to 
any new conditions. Unfortunately, that is not the case. It is a straightforward exercise to 
see that the coefficient of |1, 2 > leads to the condition 

<y{p)A{h, h)ai-h, h) = 0, (67) 

p 

where 

a{-ki, k2) = ((1 -4.g) + g (-e'^^ + 2e~'^^ + e~'^' - E^{k{))) (68) 
Similarly, at the right boundary, the coefficient of |L — 1, L > leads to 

Y,a{p)A{k^,k2)hkuk2) = (69) 
p 

where 

/3(A;i, k2) = e^('=i(^-i)+'=2(L+i)) (^^i _ 4^) ^ ^^gifca _ g-ite ^ 2e'''' - Ei{k2))) (70) 

These are new conditions which imply that the Bethe state is not an eigenstate of the 
Hamiltonian even after one finds bulk and boundary S matrices which are solutions of the 
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two body Schrodinger equations. The presence of both the momenta ki,k2 in boundary 
scattering phase shifts a, /3 in the above equations is the signature of three body interactions 
that prevent the separation of the scattering matrix into a bulk and a boundary part. The 
extra three body interaction coming from the special possibility of both the particles being 
close to the boundary thus shows us that the conventional Bethe ansatz does not quite work 
in the case of the present problem. 

Though the problem does not appear to be integrable, at least from the PABA point 
of view, we can nevertheless use the Bethe equations to construct some exact two particle 
eigenstates. This can be carried out by allowing for a linear superposition of the the Bethe 
wave functions of the kind discussed above with various values of the momenta. Instead 
of(|3Bj). we can look for a wave function of the type 

^(x, y) = -^{x, y) + T^'ix, y) + A^"(x, y), (71) 

where \E'',\E'" have exactly the same functional form as \E', but they are chosen to have 
different momenta, k[,k2, and k1,k2 respectively. The functions are also chosen to satisfy 
the scattering equations ()59|60|6H62|63|) . Furthermore, we shall choose the three pairs of 
momenta to be related by 

y: E{h) = y: E{k'^ = E E{k':) (72) 

i=l,2 i=l,2 i=l,2 

with -^(A;) given by (j47|) . We can now fix the two undertermined functions of the momenta 
A,r by requiring the state to be an eigenstate of the Hamiltonian. Clearly, the equations 
that need to be satisfied by F and A follow from ()67I69|) . They are 

Y Mh, k2)ai-k^, k2)+Tj2 (MK, k'^M-K, + A ^ {A{kl k'^)&{-kl k'^)) = (73) 
and 

Y k2)/3iku ^2) + r ^ {Aik[, k'^mk',, k'2)) + A ^ {Aikl k';)^ikl kl)) = (74) 

In the above equations the 'primed' sums refer to summations over all permutations and 
negations of momenta as we have carried out in the rest of the paper; i.e = X]p cr(p). These 
equations can indeed be solved, allowing us to construct some special eigenstates of the 
Hamiltonian. In other words, the two body problem is solvable in a sub-space of the Hilbert 
space spanned by the Bethe wavefunctions. In the presence of integrability, every Bethe state 
corresponding to any of the allowed values of momenta determined bv(|64j) would have been 
an eigenstate of the Hamiltonian. Though that is not the case, special linear combinations 
of the Bethe states do still solve the two magnon problem. It is however not clear to us if 
the above two body solutions can be utilized to solve the n > 2 magnon problem. 
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4 Epilogue: 



The test of integrability carried out in the above analysis differs significantly from the usual 
one carried out to determine whether or not particular boundary conditions maintain the 
integrability of a given integrable spin chain. One usually declares a spin chain with open 
boundary conditions to be integrable if its bulk and boundary S matrices satisfy the bound- 
ary Yang-Baxter or reflection algebra This, for example, is the notion of integrability 
used in[4 for the open so(6) spin chain. However, it has to be kept in mind that the S 
matrices satisfying the reflection algebra is a necessary and not a sufficient condition for 
integrability. Moreover, the bulk and boundary S matrices truly become matrices when 
the impurities carry some internal degrees of freedom, which is not the case in the present 
problem, where the impurity can only be of one kind. The reflection algebras too, become 
non-trivial matrix relations only when the scattering matrices describe particles with inter- 
nal degrees of freedom, and hence become non-commutative matrices. In a two component 
system like the present one, where the scattering matrices are just numbers depending on the 
momenta, the reflection equations are trivially satisfied. Thus, the analysis presented above 
also gives us a curious example of a situation where the necessary condition for integrability, 
namely the associated reflection algebra, is satisfied, while there does not appear to be a 
Bethe ansztz for the problem. 

It is also instructive to compare the present study with those concerning the Inozemtsev 
spin chain. As mentioned before, the two (and three) loop su{2) dilatation operator in the 
sector of single trace operators has been shown to be embedded in a particular integrable 
long ranged spin chain known as the Inozemtsev spin chain. The Inozemtsev chain has 
a free parameter in which the Hamiltonian can be expanded and it is part of a complete 
set of commuting charges. It has been shown in|nSlthat it is possible to fine tune this 
free parameter so that up to the three loop order, the su(2) dilatation operator and the 
Inozemtsev spin chain become identical. Indeed, this is also the easiest way to establish 
perturbative integrability up to the three loop order in the su{2) sector of the gauge theory. 
However, the picture changes when one considers non-periodic boundary conditions. To the 
best of our knowledge, the only boundary conditions other than the periodic one that are 
compatible with the integrability of the Inozemntsev and other Inozemtsev like long ranged 
spin chains such as the celebrated Haldane Shastry model are special reflecting boundary 
conditions, where the spin chain interacts with its image on the other sides of the left and 
right boundaries in very specific ways, see for example [211 12H1 I2Z1- The two loop (bulk) 
dilatation operator that we have in the present problem is the same one that arises in the 
analysis of the closed spin chains, and thus, it can be embedded in the Inozemtsev spin chain 
to this loop order. The boundary conditions in the present problem though are not of the 
reflecting type that lead to integrable 'open' versions of the Inozemtsev or Haldane Shastry 
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spin chains which can possibly be an explanation for the absence of a consistent Bethe ansatz 
for the problem. 

It is probably worth reiterating at this juncture that the two loop analysis presented in 
this paper is not the only example of the emergence of boundary conditions that lead to a 
potential violation of integrability of the dilatation operator of A/'=4 SYM. As mentioned 
earlier in the paper, the dynamical boundary conditions generated from the study of spin 
chains coupled to sub-determinant like operators in[S] produce quantum spin chains whose 
status, at least as far as integrability is concerned, is not clear. Another related scenario of 
a violation of integrability in this gauge theory arises when one studies the large limit of 
the so(6) invariant matrix model obtained as the dimensional reduction of the scalar sector 
of the super Yang-Mills theory to one (time) dimension. The quantization of the large 
limit of this matrix model was carried out recently by Klose in^Tj, and a breakdown of 
integrability of the model at higher loops was explicitly demonstrated. 

As a final observation, it might be worth comparing the present analysis to the one 
performed in the context of defect conformal field theories. M =2 defect conformal field 
theories are the conjectured gauge theory duals for the open string theory generated by 
wrapping a D5 brane on an ^^5*4 x S"^ subspace of ^^5*5 x S^. For the su{2) sector of 
such gauge theories, the open string degrees of freedom are generated by placing quark like 
fields at the ends of the spin chain, which result in open/free boundary conditions [36 } l37 | 
EHUSnHmilllllSllSlllllllSlllE]. These are different from the boundary conditions that we 
encountered in the present work, but they are similar. For such open spin chains, an all loop 
Bethe ansatz was recently proposed in[34j. Although, these are two different theories, and 
it probably does not make sense to compare their spectra, one can nevertheless compare the 
generic structure of the Bethe equation obtained in that analysis to the obstacles encountered 
in the present work. In the Bethe ansatz proposed in|^, the boundary scattering matrix did 
not have any dependence on the 't Hooft coupling of the gauge theory. Although the present 
analysis prevents us from writing down a set of Bethe equations of the present problem, we 
were able to determine the boundary scattering matrix at the two loop level, and found it 
to have a non-trivial dependence on the 't Hooft coupling. From the present analysis, it 
appears that such dependence is a generic feature of boundary scattering for (long-ranged) 
open spin chains, in the absence of reflective boundary conditions(see for example |28l I27j ). 
Thus it might be a reasonable goal to perform an explicit two loop computation for the 
5^(2) sector of M =2 defect conformal field theories, as there might be hidden surprises 
to be uncovered. At the very least, it might be possible to understand their integrability 
by studying the relation between the boundary conditions that emerge out of the defect 
conformal field theories and the reflecting boundary conditions that have been studied in 
the context of long-ranged spin chains like the Haldane-Shastry and Inozemtsev chains in 
the past|2nil2Hll2II- 
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Thus in summary, the analysis presented in this paper presents some evidence in favor 
of lack of higher loop integrability in M =4 SYM when one ventures into the study of the 
gauge theory dual of open string like operators. It is worth noting that similar obstacles in 
the direction of the implementation of Bethe ansatz methods at higher loops in the open 
string sector of this gauge theory were also arrived at recently in jlHI- The present analysis 
of integrability was carried out in a somewhat unconventional approach. It would perhaps 
be worthwhile to carry out a numerical study of the spectrum of the two loop spin chain 
presented in the paper along the lines of the involved work done in^U]. The empirical 
knowledge of the spectrum can then be used to study whether or not the degeneracies 
responsible for the one loop integrability are lifted at two loops and possibly provide a 
more direct study of integrable aspects of the problem. A numerical investigation of the 
spectrum of the gauge thoery dual of open strings coupled to non-maximal giant gravitons 
was recently carried out in|in]- In that paper some tell-tale numerical signs of integrability 
were also presented. Thus, it would be extremely interesting to understand if integrability is 
present in the problem in some way that is independent of the existence of a Bethe ansatz. 
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